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Abstract 

We present the results of the first lattice QCD study of the strong coupling gp* £>„-. 
From our simulations in the quenched approximation, we obtain qd'D-k = 18.8 ± 
2.31^0 an d 9c = 0.67± 0.08lolo6- Whereas previous theoretical studies gave different 
predictions, our result favours a large value for g c . It agrees very well with the recent 
experimental value by CLEO. g varies very little with the heavy mass and we find in 
the infinite mass limit g^ = 0.69(18). 



PACS: 12.38.Gc (Lattice QCD calculations), 13.75. Lb (Meson-meson interactions) 



1 Introduction 



Recent measurement of the full width of the charged D*-meson, T(D* + ) = 96 ± 4 ± 
22 keV Q, allowed for the experimental determination of the strong coupling of D-mesons 
to the P-wave pion gD*+D°n, namely 

g D * D7r = 17.9 ± 0.3 ± 1.9 , i.e. £=0.59 ±0.07, (1) 

where the definition of £ in |IJ is quoted below. This coupling has been extensively studied 
in the literature Q, with a variety of approaches : model independent approaches, the QCD 
sum rules, the quark models. 

Model independent approaches 0, [| have produced windows or bounds which have 
been confirmed by experiment eq. ([!]). By the way a rigorous bound g < 1 is set by the 
Adler-Weisberger sum rule j|, H. Notice that g = 1 is the naive non relativistic quark 
model result, to be lowered by relativistic corrections. 

The various QCD sum rules have been discussed with much care before and after the 
measurement in eq. ([!]) and have shown a surprising convergence towards a very low value, 
almost a factor two below eq. ([[]), see ||. In particular, in. 0, a value go'D-n = 10.5 ± 3.0 
has been quoted. No convincing explanation has been found for this discrepancy. 

Quark models, on the contrary, can accommodate large values. Good predictions have 
been produced, prior to the experimental measure || But there is a large spectrum 
of predictions 0.3-0.8, corresponding to a multiplicity of models or choice of parameters, 
therefore one can wonder whether successes are truely significant, or rather due to a happy 
choice. Without entering into details, one can answer as follows: 

1) Light-front quark models have indeed a large range of predictions 0.3-0.8, mainly 
because they use free quark Dirac spinors, which yield relativistic corrections very sensitive 
to the choice of the ill-determined light quark mass. 

2) Dirac type models, on the contrary, yield naturally large values of g > 0.6, because 
anyway, a large effective mass is generated for the light quark through the potential ||. 
These authors find g = 0.6. Too large values obtained in other calculations can be corrected 
by a quark current renormalisation factor, but one loses predictive power. 

Before claiming that the QCD based evaluations for this coupling are in conflict with 
the experimental value, it is important to compute this coupling by employing the lattice 
QCD simulations, as recently suggested in ref. [O]. An exploratory lattice calculation 



has been performed in ref. |12| but in the static limit of the heavy quark effective theory 



(HQET). Here we will directly work in QCD with relativistic propagating quarks since, 
contrary to the case of 6-physics, currently accessible lattices allow to accommodate the 
charm quark mass and therefore no heavy quark extrapolation is needed. This makes the 
lattice study of gr>*D-K rather clean. 

In this paper, we report the first calculation of this type in which we used the (improved) 
Wilson fermions. Our final result at (3 = 6.2 is 

gn*D« = 18.8 ±2.3j£J, (2) 



1 A rather exhaustive list of results for this coupling can be found in J|, To the references listed 
there, one should add also ref. p0|. 
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thus in a very good agreement with experiment. The small value predicted by the QCD 
sum rule still needs an explanation. 

The value of the coupling go*Dn provides also the access to the ^-coupling which is one 
of the main parameters of the approach based on the use of chiral perturbation theory for 
the heavy-light systems, g is related to gD*Dn through 



2m D * g c /f n 



go*DiT = { 2y/m D m D * g e /f v 



2m D g c j f % 



m 



where we use f n = 132 MeV. All the above definitions coincide up to l/m c corrections. 
Throughout this paper we will use the most symmetric definition ||. Notice that we 
assigned a subscript c, to stress that the g is not obtained with the infinitely heavy quark 
(mesons), but rather from the (not-so-heavy) charmed heavy-light mesons. By using the 
definition ||, from our result (|2]), we obtain 



9c 



0.67 ± o.o8±S:SS 



(3) 



We performed the simulations at two values of the lattice spacing (/? = 6.0 and (3 = 6.2) 
to study the systematic effects which are discussed in section 6. The simulation was firstly 
done at (3 — 6.0 and will be summarized in section 5. Our main results are given from the 
simulations at (3 = 6.2 since they have smaller 0(a)-effects. This analysis is detailed in 
sections 2 to 4. 



2 Lattice parametrization and results of the analysis 
of the two-point functions 

The main results presented in this paper and detailed in this section are obtained from 
the simulation on a 24 3 x 64 lattice with periodic boundary conditions, at /3 — 6.2. Our 
sample contains 100 independent SU(3) gauge configurations produced in the quenched 
approximation (i.e. n F = 0). The quark propagators are computed by using the following 
Wilson hopping parameters for light (q) and heavy (Q) quarks: 

K q = 0.1344 91 , 0.1348 92 , 0.1351 g3 ; 

k q = 0.1250 Ql , 0.1220 Q2 , 0.1190q 3 , (4) 

where we introduced the labels qis and Q1-3 that will be used throughout this paper. We 
have implemented the non-perturbative 0(a) improvement of the Wilson fermion action, 
by setting c S w — 1-614 [|i~4] . 



In this section we consider the standard two-point correlation functions 

cP P (t x ;p) = (£ P(0)P(x)) , C^(t x ;p) = (£ e* s V,(0)V v (x)) (5) 
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where t x > and where P = q'j^q and = q'^^q which will be used with heavy-light 
and light-light quarks. We define the constants Zp and Zy so that 

Cp> P {t x ;p) ~ Z p^j^' C^ Vi/ (t x ;p) ~ Z y 2 ^ (V - V^VvIv 2 )- (6) 

at large t x , where £p is the ground state pseudoscalar (vector) meson energy. 

From the standard study of these two-point light-light correlation functions, we ex- 
tracted the masses of pseudoscalar (amp) and vector mesons (amy), the decay constant 
(afp) and the (bare improved) light quark mass (ap) that we obtain by using the axial 
Ward identity, d^A 1 = 2pP. In table [1], we list our results for both the degenerate and 
non-degenerate combinations of our light quarks. In the computation of the pseudoscalar 





amp 


amy 


ap 


af* 


K qi 




0.306(1) 


0.409(3) 


0.0413(4) 


0.068(2) 


K qi 


~ K Q2 


0.284(2) 


0.393(4) 


0.0354(3) 


0.065(2) 


K qi 


— K qs 


0.266(2) 


0.381(4) 


0.0310(3) 


0.063(2) 


K q2 


— K q2 


0.259(2) 


0.377(5) 


0.0296(3) 


0.062(2) 


K q2 


— K q3 


0.240(2) 


0.364(5) 


0.0252(3) 


0.060(2) 


K q3 


— K qs 


0.219(2) 


0.350(6) 


0.0208(2) 


0.058(2) 



Table 1: Light meson masses, bare quark masses and (renormalized) pseudoscalar decay 
constants. The time intervals chosen for the fits are: P : t G [10,30], V : t £ [11,25], and 
p,fp:te [12,29]. 



decay constant and of the bare quark mass we improved the axial current at 0(a), i.e. 

A[(x) = A,(x) + c A (g 2 )d ll P(x) , (7) 



where ca = —0.038(4), as determined non-perturbatively at (5 = 6.2 in refs. |14|, [H| [16 



In the computation of the renormalized decay constant, the 0(ap) effects are eliminated 
by redefining 



Zi{9l) = zf(gl)(l + b A (gl)ap) , 



where the non-perturbatively estimated constants are Z A °^ = 0.81(1) |L5], [L7|, jlSf , and 
b A = 1.19(6) U 
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By using the lattice plane method illustrated in fig. ffl, we get af n = 0.0488(24), 
which after comparison with the physical f n = 0.132 GeV, leads to the following value of 
the inverse lattice spacing: 

a-\^) = 2.71(12) GeV. (9) 

That value is consistent with the one obtained by using the p- meson mass (a~ 1 (m p ) = 
2.62(9) GeV) and/or the K*-meson (a- l {m K *) = 2.67(8) GeV). 

As for the heavy-light systems, in what follows, we will need the pseudoscalar and 
vector meson masses, in addition to the constants Z P , Z v ©. We obtain those quantities 
by fitting our lattice results for the correlators (^j) with the mesons at rest (p = 0), to the 
forms given in eq.(^). The results are presented in table 0. 





amp amy Zp Zy 


K Qi ~ K qi 

K Ql ~ K q2 
K Ql ~ K qz 


0.692(2) 0.732(2) 0.0187(9) 0.0068(4) 
0.678(2) 0.719(3) 0.0176(9) 0.0062(5) 
0.669(2) 0.710(4) 0.0171(10) 0.0059(6) 


K Q2 ~ K qi 
K Q2 ~ K q2 

K Q2 ~ K qz 


0.788(2) 0.822(2) 0.0207(10) 0.0080(5) 
0.775(2) 0.809(3) 0.0195(11) 0.0073(6) 
0.766(3) 0.800(4) 0.0189(12) 0.0070(7) 


K Qz ~ K q\ 
K Qs ~ K q2 
K Qz ~ K qz 


0.878(2) 0.907(2) 0.0222(11) 0.0091(6) 
0.866(2) 0.894(3) 0.0208(12) 0.0083(7) 
0.857(3) 0.885(4) 0.0203(13) 0.0079(8) 



Table 2: Meson masses and Z's in lattice units extracted from our lattice data. Indices in the 
hopping parameters (qis and Q1-3) are specihed in eq. (0). The time-intervals used for the fits 
are P : t G [15,30], V : t G [17,28]. 



3 Computation of the three-point functions: Basics 
and Results 

In this section we define the quantities that we need to compute and explain the strategy 
that will allow us extracting the coupling gyp^, where P and V stand for the heavy-light 
vector and pseudoscalar mesons respectively. We will then present our results obtained for 
the heavy and light quark masses that are directly accessible from our lattice. 
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3.1 Theoretical basis 



We start from the computation of the transition matrix element between the heavy-light 
vector meson (V) and the heavy-light pseudoscalar (P), mediated by the axial light-light 
current = tpf^q (see fig. |2|). It is parametrized as 



(P(j/)\A»\V(p,\)) = 2m v A (q 2 ) € -^-q» + (m P + m v )A 1 (q 2 ) 



+A 2 (q 2 ) 



e -q 



ITLp + ITLy 



p v _|_ j/i* _ m v - 



2 ™2 

2 <r 



(10) 



where q = p — p' 0. The matrix element of the divergence q^A^ is dominated by the pion 
pole for q 2 in the vicinity of m 2 : 



{P{p')\q^\V(pA))=9vP, 



e A (p) ■ g 

9 9 
m4 — ? 



x /,mj + 



From eq. (11), we find for q = : 



9vp-k 



2m v A (0) 



(12) 



With our settings, the lattice transfer at q = happens to be close to q 2 = (see table §). 
However, the lattice simulations at q = can only give the form factor A\. The other ones, 
y4 0j 2, can be computed at q 7^ which in our case is no longer close to q 2 = 0, as can be 
seen in table [5| Since has the pion pole, it varies very fast in the vicinity of q 2 = 
and thus cannot be directly extrapolated. To overcome this difficulty, we express A (0) in 
terms of y4 12 (0) which do not have a pion pole. This relation can be obtained by using 
the fact that in eq. (|T0"D, the axial current cannot have a singularity at q 2 = 0. The three 
residues of l/o 2 -terms must cancel, this leads to: 



9vPn — -r~ 



(m v + m P )A 1 (0) + (my - m P )A 2 (0) 



(13) 



In this relation, one can see that the A\ contribution is dominant. This dominant con- 
tribution is directly obtained from lattices at q — 0, as already mentioned. A2 has to be 
extrapolated from q 7^ to q 2 = 0. Its nearest pole is the a\ meson mass. The error 
generated with this extrapolation will be discussed in the next subsection. This error is 
anyhow harmless since it applies only to the subdominant contribution (< 5%) to gyp-ir- 

In the heavy quark limit (itlq —>■ 00), in which my = mp, only the form factor Ax(0) 
contributes. In that limit one recovers the formula used in ref. |12| to compute this coupling. 
For heavy quarks close to charm, the corrections in powers of the inverse heavy quark 
mass are expected to be sizable, which is why we decided to compute gvPn with the 

2 This definition of form factors is equivalent to the usual one : the matrix element (V P) has a 
"-" sign in front of A2 . 



6 



relativistic (propagating) heavy quark. Notice that isospin symmetry relates various charge 
combinations of the VPtt couplings: 

gvPn = gv+p°n+ = -V2g v + P + n o = -g v op+ n - ■ (14) 
For simplicity, we define 

2 m v + m P 2 2 m v -m P 2 
Gi{q ) = A^q ) , G 2 {q ) = - A 2 (q ) , (15) 

Jn Jn 



and rewrite eq. (0) as 



As already mentioned G\(0) is the dominant contribution to gvPn, the G2/G1 being a few 
percent correction to it. 

Analogously, the expression for gg at a given heavy quark Q mass is 



where we note, according to 



yQ 2y/m v m P v ; 



3.2 Computation of the three point functions 



To access the matrix element (|10f) from the lattice, we compute the following three-point 
functions 



C^(0,q,t x ,0,t y )= (J2 V fl (0)A u (x)P(y)e- i ^) 



(19) 

0<t X <ty 



where the pseudoscalar meson is inserted at rest {p ' = (0,0,0)) at a fixed time chosen to 
be t y = 31a, while the current operator has a momentum q G {(0,0,0), (1,0,0)} in units 
of the elementary momentum (2ir/La ~ 0.7 GeV). The vector meson interpolating field is 
at the origin (0, 0) of the lattice. 

When both mesons are at rest, the only useful ratio is (t = t x ) 

Ri(t) = { } ( l p , (20) 

which develops a plateau for t G [12,17]. At that plateau we extract the matrix ele- 
ment (PHI), i- e - the value of the form factor A 1 (q 2 ). To access the ratio A 2 /Ai, we study 
the ratios with the momentum injection q = (1, 0, 0) x 2n/L, namely 

^ ^ ~ ^(2) 77 ^(2) 



Cy 2 v 2 (^ ^)Gpp(t y — t) 
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Cy 2 V 2 (t'> 0)Cpp(t y — t) 



r (3) (t-a)Z 1/2 Z 1/2 
m = ^^% Zp ■ (21) 

After inspecting the ratios ([H]) for all the combinations of heavy and light quarks, we choose 
to fit them to a plateau at t 6 [14, 17]. We illustrate in fig. § the signals for all the four 
ratios for a given kq, tz q . We notice that the ratio R2(t) has a hardly observable plateau, 
the value of which is compatible with zero. As just stated the ratio R\(t) contributes more 
than 95% to gyp p and its plateau is very good. 

Denoting r; the average value of Ri on the plateaus, the form factors A% 2 are then given 
by: 

A 1 (q = 6) = ^ ' A 1 (q=^(l,0,0))= ^ , (22) 

my + mp \ L ) my + mp 



Ax 



(mp + m v ] 



(q 2 - Ey(Ey-mp)) + 



2mpq 2 

where Ey is the energy of the vector meson. 



m v (E v -mp)r3 m v q 1 r 2 

Ey 7 4 +l ^y~n 



,(23) 



Our numerical simulations indicate that the ratio A2/A1 is positive and of the order of 
1. It results from eq. ( p!5|) that G2/G1 ~ (my — mp)/(my + mp) which leads to a small 
and positive correction to Gi(0) in (|16|) . 

In the left part of table |3|, we present our results for the form factor Ai(q 2 ) (Gi(q 2 )) for 
all the quark combinations and with both mesons at rest, q 2 = (my — mp) 2 . In this case, 
q 2 E (0.04, 1.16) • 10~ 2 GeV 2 . They are obviously very close to zero and it is reasonable to 
assume that, at such small q 2 , G\(q 2 ) ~ Gi(0). 

The results for the ratios of the form factors AijA\ (or, equivalently, G2/G1) with 
q = (1,0,0) x 2n/L are listed in the right part of table |3[ They are obtained from the 
ratios r 2 _4 (see eq. (p3|) ). As expected, the ratios of the form factors G2/G1 are positive and 
very small (they never exceed 5%). Now this ratio has to be extrapolated to q 2 = 0. The 
Ai }2 form factors may change significantly but without changing sign since the nearest pole, 
the ai-meson, lies at larger q 2 (m ai ~ 1.2 GeV |B|) Q. On the other hand, the ratio A 2 /Ai 
is expected to be rather constant because A li2 have the same pole factor 1/(1 — q 2 /m 2 ai ) 
which cancels out in the ratio We thus extrapolate A 2 /Ai by keeping this ratio constant. 



3 This is confirmed by the ratio r 4 /n = A 1 (q= (1,0,0) x 27r/L)/A 1 (q= 0) - 0.6. 

4 1/(1 - q 2 /m 2 ai ) varies of about 30% - 40% since for q = (1,0, 0) x 2tt/L, q 2 E (-0.48, -0.44) GeV 2 . 
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q 


= (0,0,0) 




q = 


(1,0,0) x 2n/L 




q 2 x 10~ 2 GeV 2 


Ax{q 2 ) 


Gi(g 2 ) 


q 2 x 10" 


2 GeV 2 Aij Ax 


Gil G\ 


K Qi 


— Kqi 




1.17(17) 


0.71(4) 


14.9 ±0.9 


-44.8 


v3-9; 


0.85(14) 


0.024(4) 


K Qi 


— K qi 


~~ K Q2 


0.52(12) 


0.68(6) 


14.8 ± 1.2 


-46.3 


^4.0; 


0.74(12) 


0.014(2) 


K Qi 


— Kqi 


— K q3 


0.23(10) 


0.68(6) 


15.1 ± 1.3 


-47.1 


(4.1; 


0.70(19) 


0.009(2) 


K Qi 






1.20(22) 


0.69(5) 


15.5 ± 1.2 


-44.7 


(3.9; 


0.93(22) 


0.027(6) 


K Qi 


— K q2 


— K q3 


0.73(21) 


0.67(7) 


15.5 ± 1.5 


-45.7 


(3.9; 


0.85(25) 


0.019(5) 


K Qi 


— n q3 


- Kq 3 


1.24(30) 


0.68(5) 


16.0 ± 1.5 


-44.5 


v3-8; 


1.13(35) 


0.034(9) 


K Q 2 


— Kqi 


— K qi 


0.82(12) 


0.71(4) 


16.9 ± 1.1 


-46.1 


(4.1; 


0.89(14) 


0.018(3) 






— K q2 


0.31(9) 


0.68(6) 


16.7 ± 1.4 


-47.3 


(4.2; 


0.77(14) 


0.010(2) 


K Q 2 


— Si 


— K q . 3 


0.11(7) 


0.68(6) 


17.1 ± 1.5 


-48.0 


(4.2; 


0.73(22) 


0.006(2) 




— K q2 


— Kq 2 


0.84(17) 


0.69(5) 


17.5 ± 1.5 


-46.0 


v 4 -o; 


1.19(22) 


0.025(4) 


K Q 2 


— K q2 


— K q3 


0.47(16) 


0.68(7) 


17.7 ± 1.7 


-46.9 


(4.1; 


0.92(29) 


0.015(4) 


K Q 2 


— K q3 


— K ga 


0.85(24) 


0.68(6) 


18.3 ± 1.8 


-45.9 


^3.9; 


1.27(38) 


0.028(7) 


K Qs 


- K qi 


~ «gi 


0.60(9) 


0.71(4) 


18.9 ± 1.3 


-46.9 




0.93(14) 


0.015(2) 


K Qs 


— K qi 


— K q2 


0.19(6) 


0.68(7) 


18.5 ± 1.7 


-48.0 


^4.2; 


0.80(15) 


0.007(1) 


K Q 3 


— K qi 


— K q , 3 


0.04(5) 


0.68(7) 


19.1 ± 1.8 


-48.7 


^4.3; 


0.78(25) 


0.003(1) 


K Qs 


— K q , 2 


— K q2 


0.61(14) 


0.70(5) 


19.6 ± 1.7 


-46.9 


(4.1; 


1.08(23) 


0.018(4) 


K Q 3 


— K q2 


— K q3 


0.62(20) 


0.68(7) 


19.9 ±2.1 


-47.7 


^4.2; 


1.00(32) 


0.012(3) 


K Q3 


— K q3 


— K q3 


0.25(8) 


0.68(6) 


20.5 ±2.2 


-46.9 


(4.1; 


1.43(41) 


0.024(6) 



Table 3: Ax, Gx,Ai/Ax,Gi/Gx at (3 = 6.2 for several values of q 2 in the physical units (GeV 2 ) 
as obtained by using a -1 = 2.71(12) GeV (from eq. (||)). 
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3.3 Comment about the effect of improvement on the form fac- 
tors 

As a side remark, we comment on the effect of the improvement of the bare axial current 
(see eq. (|7])) on the form factors Ai 20 (q 2 ). The divergence of the current in the matrix 
element ([LO]) leads to 

(P(p')\qjsq\V(p,\)) = 2m v { -^- A (q 2 ) , (24) 

where we used eq. ( PH ) and the axial Ward identity, d^A 1 ^ = 2p q P. Therefore the improve- 
ment will only affect the form factor A (q 2 ), but not the other two (Ai t2 (q 2 )). In other 
words, 

A li2 (q 2 ) - A^q 2 ) = A li2 (q 2 ) , 
Ao(q 2 ) - Kiq 2 ) = (l - c A ^A (q 2 ) . (25) 

As our method to compute gD*Dn relies on the computation of the form factors A\ and A 2 
only, our result does not depend at all on the improvement of the bare axial current. 



4 Chiral extrapolations and the heavy quark interpo- 
lation 

Now it is a simple matter to combine our results from table ||] in the way indicated in 
eq. (|16|), and to compute the values of the couplings gvp P f\ Those numbers should now 
be extrapolated to the to the physical pion mass, to get gvPn for each of our heavy quarks. 
Finally this is to be followed by interpolation in the heavy meson masses to the ones 
corresponding to D and D*. 



4.1 Chiral extrapolations 

Our light pseudoscalar mesons ("pions") are rather heavy and we need to extrapolate 
to the physical pion mass [], m n = 0.14 GeV. From the lattice planes method pl]| , the 



physical pion mass in lattice units is am^ = 0.053(2). To obtain gvpn, we can fit the 
5VPp-dependence on the quark mass with three formulae: 

- Linear extrapolation: by fitting our data to 

gvPp = a>o + ai(am p ) 2 , (26) 

where the values of (am p ) are listed in the first column of table |l|. This fit allows 
to fix the parameters a 0) i and thus to extract gyp-K by either simply reading off a 
(chiral limit), or by combining a 0j i with am n . We will give results by opting for the 
latter choice. 



5 The lower subscript "p" labels the light pseudoscalar meson. 

6 In physical units, our light pseudoscalar mesons are in the range m p £ (0.5,0.8) GeV. 



10 



- Quadratic extrapolation: we also attempt the quadratic fit 



9vp p = b + bi(am p ) 2 + b 2 (am P Y 



(27) 



"Chiral Log" extrapolation: One can also fit to a form motivated by the chiral per- 
turbative expansion for the heavy-light systems (for a review, see ref. In par- 

ticular, the coupling of the pion to the heavy meson doublet receives logarithmic 
corrections [E2l : 



9vp p = c + Ci(am p ) 2 + c 2 (am p ) 2 \og((am p ) 2 ) . 



(2f 



The results of all three fits are presented in table |] and illustrated in fig. f|. The 
quadratic and logarithmic extrapolations are rather unstable i.e. significantly dependent 
on small changes in the analysis procedure, because we extrapolate from too heavy quark 
masses. We present them mainly as an estimate of the systematic error. 



amp 



am v Gi(0)( ,in -) G^O)^-) G^O)^- 



(lin.) 

9vptt 



Qi : 0.1250 
Q 2 : 0.1220 
Qi : 0.1190 



0.645(3) 0.688(6) 16.8 ± 1.9 15.6 ±2.4 15.0 ± 3.0 17.7 ±2.2 
0.744(4) 0.779(6) 19.3 ± 2.4 18.0 ± 3.3 17.2 ± 4.2 20.1 ±2.7 
0.836(4) 0.865(7) 21.7 ±3.0 20.3 ± 4.3 19.3 ± 5.6 22.6 ± 3.3 



Table 4: For each heavy quark directly simulated on the lattice we show the values of the heavy- 
light meson masses for which the light quark is (linearly) extrapolated to the u/d-quark mass. 
We also list the values of the corresponding (?i(0) (i.e. the gvPn coupling without the G 2 /Gi 



correction) by using eqs. (2t, |2?J 28), and in the last column, the gvPn linearly extrapolated 
including the G 2 /G\ correction (see eq. (|l6p). 





^(O.lin.) ^(lin.) 

9q 9q 


Q x : 0.1250 
Q 2 : 0.1220 
Qi : 0.1190 


0.636(65) 0.669(72) 
0.641(76) 0.668(79) 
0.648(86) 0.673(88) 



Table 5: For each heavy quark we show the values ofg^ for which the light quark is linearly 
extrapolated to the u/d-quark mass and the gQ linearly extrapolated including the G 2 (Q) j 'Gi(O) 
correction (see eq.dT^)). 



11 



4.2 Chiral extrapolations with two different light quarks 

We have computed the three-point Green functions with two different light quark masses 
1711,1712 for the quarks (71,(72 shown in fig ||]. The form factors extracted from the latter 
depend, in general, both on 1712 + vii\ and m>2 — mi what makes the chiral extrapolation 
rather tricky. Luckily, however, it can be shown [] that, in the heavy quark limit, the 
dominant contribution to gypn, Gi(0), only depends on m 2 + m\. In this limit, 
we can use data with (rrii 7^ m 2 ) to perform our chiral extrapolation as a function of 
ni2 + mi (or similarly, as a function of (m p ) 2 ). We have compared the results of the 
linear extrapolation using this method with the one which uses Green functions with only 
rrii — m 2- The extrapolated results agree within 2% while the statistical error is typically 
~ 15%. Considering six (m^mj) couples allows to estimate the systematic error due to 
the chiral extrapolation by comparing the linear, quadratic and logarithmic extrapolations 
as explained in the preceding subsection. As can be seen in table f|, Gi(0) was computed 
in this way. 

The quark mass dependence of the corrective term G2(0)/Gi(0) in eqs. ( P~0| , |TT[ ) is dom- 
inated by the factor my — mp in eq. (|T5|). In the infinite mass limit, it is known that 
my — mp oc m 2 — mi + O{\jmo). The corrective term G2(0)/Gi(0) is then dominantly 
proportional to m 2 — mi and the contributions with mi 7^ m2 are useless for the chiral 
extrapolation. Therefore, we have only considered the case mi = m-2 when including the 
correction G 2 (0)/Gi(0). 



4.3 Interpolation to the charm sector 

Finally, we need to reach the mass of the charm quark. To that end we will use the values 
of the spin-averaged masses of the heavy-light mesons 

_ 3m v + m P 

m H = . (29) 

When converted to the physical units by means of a~ 1 (f 7T ) (given in eq. (Q)), we have 

m H e { 1.83(9), 2.08(10), 2.32(11) } GeV . (30) 

and mp, = 1.974 GeV is within this range. The coupling 7)q that we already mentioned in 
the introduction : 



2 v /m P m y _ 

QvPir = 7 9Q ■ (31) 

Jit 

is the proper parameter to be used for an interpolation in the heavy quark mass. The 
heavy quark symmetry suggests us to fit our data to the following forms 

9vpt7 a 2 ^ , 62 , qo n 
— — = ai + ^— and go = bi + =r - . (32) 
mn mpi mn 



7 Using the heavy quark symmetry and the hermiticity of A^, one derives that 
< V{m 2 )\A tl \P{m 1 )>=< P{m 2 )\A tl \V{m 1 ) >=< V(mi)\A^\P{m 2 ) >■ 

8 We have indeed checked that the slope of G2(0)/Gi(0) as a function of 7712 — mi is four times larger 
than the one as a function of mi + m 2 . 
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Using these forms we now interpolate to the charm region each of our three sets of chirally 
extrapolated data from the previous subsection. The results are presented in table || 



X-extrap. 


9D*Dit 9c 9oo 


linear 
quadratic 
X-log 


18.82 ±2.34 0.669(75) 0.69(18) 
17.11 ±3.10 0.574(88) 0.62(27) 
16.39 ±4.02 0.551(118) 0.57(37) 



Table 6: Results of the linear interpolation of the form (p^j to D-D* mesons at (3 = 6.2, for each 
of the three chiral extrapolations discussed in the text. Note that only the linearly extrapolated 
result incorporates the G2/G1 correction. The other two are uncorrected. 

Illustration of that interpolation for the case of the linear chiral extrapolation is provided 
in fig. H Note that our results indicate that the slope in 1/mg for the coupling gQ is 
small and negative. Assuming that the linear dependence in I/tuh holds all the way to 
1/rriH 0, we obtain that ^ is not more than 15 % larger than g c (numerical results for 
(Too are also given in table 

5 Results at /3 = 6.0 

In this section we briefly summarize the results obtained at (3 = 6.0. The analysis fol- 
lows the lines presented in the previous sections. We have run in a 16 3 x 64 volume 
over 100 configurations with the following set of light quarks Wilson hopping parameters: 
n q = 0.1339, 0.1342, 0.1344, 0.1346. We will only show results for one heavy quark, the 
closest to the physical charm, kq = 0.1190 (mp = 1.77(11) GeV), for several values of the 
light quark mass (see table 0). For (3 = 6.0, the ratio G 2 (0) / Gi(0) has about 100% error, 
and we prefer to give the uncorrected result (see eq. (0)), in table [5]. 

It is worth noticing that the results at /3 = 6.0 agree within errors with those at 
(3 = 6.2. This is illustrated in fig. || where we compare the case j3 — 6.0, kq = 0.1190 to 
(3 = 6.2, kq = 0.1250, which corresponds to approximately the same physical mass of the 
heavy-light meson (mp = 1.77(11) GeV and 1.75(9) GeV respectively). Notice that the 
numerical results at (3 = 6.0 are more noisy than those at /3 — 6.2. 

For this reason and also because the lattice spacing at /3 — 6.0 implies strong limitations 
on heavy quark masses, we did not attempt an extrapolation to the charm and even less to 
the infinite mass limit. Notice nevertheless that the charm mass region is almost reached 
with kq = 0.1190 in our setup. 

6 Physical results and discussion of errors 
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0.1339 


0.1342 


0.1344 


0.1346 




Gx(0) 
9q 


19.9 ±2.4 
0.80(7) 


20.8 ±3.0 
0.80(9) 


22.7 ±4.3 
0.83(12) 


19.3 ±4.0 
0.78(13) 


21.4 ±6.0 
0.80(19) 



Table 7: Values of Gi(0) and gk at (3 = 6.0 for one fixed heavy quark: kq = 0.1190 and 
several light quarks. The extrapolation to the physical pion mass is presented in the last column. 

6.1 Systematic uncertainties 

• Discretisation errors- 1: In our study we implemented the full 0(a) improvement of 
the Wilson QCD action and the axial current. As we discussed in the text, the 
improvement of the bare axial current does not influence the value of our gD*Dn- As 
for the renormalization constant, we used the non-perturbatively determined value, 
including the coefficient b^, which ensures the elimination of the artifacts of 0(ap). 

• Discretisation errors-IT. Our calculation has been made at (3 = 6.2. With intention 
to study the 0(a) effects, we have also performed the simulation at (3 — 6.0 which 
has been summarized in the preceding section. The good agreement illustrated in 
fig. |5] shows a small discretization error. One might wonder if the small decrease of 
g from (3 = 6.0 to (3 = 6.2 is the sign of a systematic finite a effect. This would then 
point toward a continuum limit lower than the values quoted here. This difference 
might also simply be a statistical one since it is smaller than one standard deviation. 
With only two values of the lattice spacing it is not possible to try a systematic study 
of the continuum limit and thus to discriminate between these hypotheses. Further 
studies at different values of (3 are badly needed. From the experience about similar 
quantities one might hope that the continuum limit will not be too different from the 
result at f3 = 6.2. 

• Discretisation errors-Ill: When interpolating to the charm quark, i.e. to the D- 
meson (|32|), we used the value of the lattice spacing as obtained from the pion decay 
constant. If the lattice spacing is fixed by the p-meson mass, the value of the go*D-K 
remains practically unchanged. This is not surprising since the slope of gvPir/ m H in 
1 /rriH is very small so that the slight change of the position of the 1/fho does not 
make any visible impact on the final result (it increases by ~ 1%). 

• Chiral extrapolations: This source of uncertainty actually dominates our systematic 
error bars. This shows up also in the difference between linear, quadratic and log- 
arithmic fits. Indeed, at (3 — 6.0 the additional lightest quark (k = 0.1346) has 
been added in order to reduce this error. We decide to take the result of the lin- 
ear extrapolation as our central value since our data follow rather clearly that form. 
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The difference between that central value (of the linear extrapolations) and the ones 
obtained through the quadratic and the "log" fits will be included in the systematic 
uncertainty. 

Finite volume effects: We also studied the finite volume effect by performing two 
parallel simulations at (3 — 6.0, with lattice of size 16 3 x 64 and 24 3 x 64. To illustrate 
the net effect, we plot in fig. [7| the ratio of a three-point correlation functions ( |19"D 
as obtained from the simulations with two lattice volumes and for both heavy-light 
mesons being at rest. Within our statistics, we do not see any evidence for the 
presence of finite lattice volume effects. To be conservative, however, we will take 
into account the observation that the central values are in the interval 



cf (t ; 24 3 ) 

and thus will include 6 % in the systematic error. 



0-95 < t 3) \'„ J < 1-07, (33) 



Quenching effects: When discussing the chiral extrapolations we also considered the 
effect of using the leading chiral log behavior. It is important to stress that such a 
behavior is only valid for the full (unquenched) QCD. In the quenched approxima- 
tion, however, one encounters the so-called quenched logs which are not of the form 
mplog(mp) as in full QCD but rather divergent, of the form m^Xogimfy . Here m 
stands for the mass of the 7/'-meson, which in the quenched theory does not decouple 
from the octet of light pseudoscalar mesons. One could thus envisage a fit to the 



form p3 | 



9vp p = Co + ci(am p ) log(am p ) + q '' log(am p ) + c 2 (am p ) +... (34) 

If we had been in the region of very small masses we should have used this form 
of the fit to determine the coefficient c ^ uench ^ anc i then correct for it when getting 
the quenched physical results. In our case, however, the meson masses we were able 
to simulate directly are in the region in which the dependence of gyp v on (am p ) 2 is 
linear as can be seen in fig. 4. 

6.2 Conclusion 

In this paper, we made the first lattice QCD study of the strong coupling of the pion to the 
spin-doublet of D-mesons, Qd*d-k- Our results, obtained in the quenched approximation, 
are in a very good agreement with the recent experimental measurement |l| . Our numbers 
are much larger than the predictions made by using various QCD sum rule techniques. It 
is also larger than the previous lattice estimate which has been made in the static limit 
of the heavy quark effective theory Jl2| and on a very coarse lattice. The reason for that 



disagreement remains to be understood. It should be noted that our data suggest that the 
dependence on the heavy quark mass is very weak and that the static value g^ > 7j c . A 
careful extrapolation to the continuum limit is still needed but, inspired from the experience 
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with similar quantities, one might hope that the continuum limit will not be too different 
from the result at (3 = 6.2. 

On the basis of our results and from the discussion of the systematic uncertainties which 
we combine in the quadratic sum, we conclude that 

gD*D* = 18.8 ± 2.3±|;J , and g c = 0.67 ± 0.08±g;| , (35) 

which is the result that we quoted at the beginning of this paper (^].|3]). 

The coupling 7j varies very little with the heavy mass and we find in the infinite mass 
limit 

(Too = 0.69 ±0.18. (36) 

A further improvement of our results includes the increase of the statistical sample, 
computations with static heavy quarks, computations at larger (3 (such as /3 = 6.4), and 
most importantly the attempt of an unquenched study of this coupling. 
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Figure 1: We illustrate the so-called method of lattice planes by which we get the value of 
a fK/ir : I n the upper plane, the points in which the fit (solid) line to our data crosses the dashed 
lines, corresponding to amy = C K ^\J (amp) 2 with Ck = (mx* /mK)phys an d C n = (nip/m^phys 
respectively, determine the values of [amj^j^) 2 , denoted by diamonds. These values are then used 
to fix af K / n (filled squares) in the lower plane, where we fit the pseudoscalar decay constants as 
a + P(amp) 2 . 
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Figure 2: The graph of the three-point function that we compute in this work. 
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Figure 3: Signals for the ratios Ri-i(t) (real part of Ri,R$,R4 and imaginary part of R2), 
defined in eqs. ( pc| |2l| ), as computed on our lattice. Illustration is provided for kq = 0.1220 and 
K q = 0.1348. 
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Figure 4: Chiral extrapolation of gyp v without the GijG\ corrections for a fixed heavy quark 
at (3 = 6.2. The curves obtained from the fit to eqs. 

dH3>H) are displayed. The filled square 
point with its error bars corresponds to the linear extrapolation. The complete list of result can 
be found in table |4j. 
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Figure 5: Fit of our data (empty symbols) for gypn &nd gQ, at (3 = 6.2, to the forms (|3^), The 
results of interpolation to the D-meson is denoted by the filled symbols. 
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Figure 6: Chiral extrapolation of g^ for a fixed heavy quark kq = 0.1190 at (3 = 6.0 and 
kq = 0.1250 at /3 = 6.2 corresponding approximately to the same heavy meson mass. 
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Figure 7: Ratio of the three point functions (19) as obtained for the same meson masses at 
j3 = 6.0 but on different volumes 16 3 x 64 and 24 3 x 64. Illustrated is the case of the heavy quark 
corresponding to kq = 0.1220 and the light ones to n q = 0.1344. 
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